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Abstract 

A recent interpretation of the caloric curve based on the expan- 
sion of the abraded spectator nucleus is re-analysed in the framework 
of the Time-Dependent Hartree-Fock (TDHF) evolution. It is shown 
that the TDHF dynamics is more complex than a single monopolar 
collective motion at moderate energy. The inclusion of other impor- 
tant collective degrees of freedom may lead to the dynamical creation 
of hollow structure. Then, low density regions could be locally reached 
after a long time by the creation of these exotic density profiles. In 
particular the systematic of the minimum density reached during the 
expansion (the so-called turning points) appears to be different. 
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1 Introduction 



During the past decade, many works both in experimental and theoretical 
nuclear physics, have been devoted to the search of the liquid-gas phase 
transition in Heavy-Ion collisions at intermediate energies. One of the most 
striking results is the observation of the so-called caloric curve 0. Indeed, 
using an isotopic thermometer in conjunction with a measure of the total 
excitation energy, the authors of ref. |3J have observed that, over a wide 
range of excitation energy, the temperature of the abraded source seems con- 
stant while for a system in an unique phase, one would expect a monotonic 
increase of the temperature with excitation energy. This curve is similar to 
those known for liquid-gas phase transition at constant pressure. Indeed, in 
a coexistence region of a single fluid system, it exists a univocal relation be- 
tween all intensive thermodynamical quantities. From the experimental point 
of view this curve is actually subject of debates since neither the tempera- 
ture nor the excitation energy are well defined quantities for small objects in 
rapid evolution such as hot nuclei @. 

Despite all these difficulties, this observation might be a major step for- 
ward since it could be a direct evidence of the liquid-gas phase transition 
in nuclear matter [El. Recently, an interpretation of this curve has been 
proposed by Papp and N6renbergj4|. Considering the expansion of excited 
projectile like fragments in a collective model, they predict that the temper- 
atures at maximum dilution follow approximately the experimental caloric 
curve. Indeed, the temperature remains almost constant while more and 
more excited systems are considered because the densities at the turning 
point of the monopole oscillation rapidly decrease as a the initial excitation 
increase. As stressed in a recent review article 0, this interpretation is of 
great importance since, it gives an intuitive connection between the creation 
of compound nucleus, radial flow, multifragmentation and the final results of 
Heavy-Ion collisions. 

In this paper we would like to re-investigate the same scenario within a 
microscopic approach: the time-dependent Hartree-Fock (TDHF) dynamics. 
This is not the first attempt in this direction, since TDHF|| [5], H as well 



as BUU||, [T(J have been used in the past to study expanding nuclear sys- 
tems. However, in this article, we will concentrate on the possible differences 
between a full quantum treatment and a macroscopic collective modelQ. 

In the next chapter, we will recall the general framework of the work done 
in ref. [Q. In part 3, we describe features of the microscopic model we have 
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used. Two forces of Skyrme type are used, one leads to a Soft Equation of 
State (EOS) in the infinite nuclear matter, the other one to a stiff EOS. These 
EOS are key issues when discussing liquid-gas phase transition and in par- 
ticular mechanical instabilities occurring inside the coexistence zone, namely 
spinodale instabilities. In this article, we concentrate on the dynamics of 
expanding systems as predicted by the time-dependent Hartree-Fock theory. 
In order to discuss the observed results we have defined an equivalent of the 
EOS for finite systems. We have pointed out many differences between the 
energy of finite systems as a function of the dilution and the infinite nuclear 
matter case, this could be in part understood in terms of surface effect and 
Coulomb interaction. 

We have paid a particular attention in initial conditions, in order to be 
in same framework as in ref. The authors of ref. [Q follows the idea 
of a dominant monopolar collective mode. This idea was used many times 
in Heavy-Ion Physics at intermediate energy in order to have quantitative 



informations on how an excited expanding uniform sphere could break [1 1 



Our calculation do not make this assumption, since we do not impose the 
domination of one collective mode on the others |12|. In this more general 



framework, the expansion dynamics appears more complex than the one de- 
scribed in ref. 0]. In particular, for moderate excitation energies (T < 
5 MeV) we show that hollow structures can be formed during the TDHF 
dynamics after several monopolar oscillations. These exotic shapes are not 
included in breathing mode pictures |IT] and could change our physical under- 



standing of Heavy- Ion collisions at intermediate energy. In a recent work[13 



we have shown, that other collective degrees of freedom should be included 
in order to describe the TDHF expansion. In the following, we will generalize 
the simple breathing mode picture in order to take into account more col- 
lective variables. In particular, we illustrate the inclusion of more accessible 
density profile by defining a generalized EOS that include hollow density pro- 
files. On the other hand, considering highly excited sources (T > 5 MeV) we 
observe that the quantum dynamics differs from collective models because of 
the importance of wave propagation and of thermal mixing on the dynamics: 
barrier transmision, large amplitude motion, lost of collectivity, ... 



2 One model for the interpretation of the 
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caloric curve 



Let us first recall the scenario followed by Papp ad Norenberg in order to un- 
derstand the ALADIN results concerning the fragmentation of a gold projec- 
tile after its interaction with a gold target at 600 MeV/A incident energy Jl]. [4| . 
Since the experiment is observing the decay of the projectile-like fragments 
in peripheral heavy ion reactions, the initial conditions are provided by the 



abrasion-ablation model |L4||: after collision, the projectile of mass A has 
lost AAi nucleons; the initial excitation E* of the nucleus is assumed to be 
related to AAi by 

E* = 7 AA (1) 

In the following simulations, AAi varies from 8 to 108 nucleons and 7 = 
13.3MeV. Temperatures are adjusted in order to get the excitation energies 
E* (H]) for a spherical nucleus of mass A- AAi. The initial density is taken to 
be lower than the normal density, = 0.8 po because the abrasion-ablation 
dynamics is assumed to induce a small dilution through the ejection of few 
nucleons 0]. 

Then, authors of ref.Q] follow the expansion of the system as an isen- 
tropic self-similar monopole motion: the breathing mode. In this picture, 
the internal degrees of freedom are assumed to be equilibrated like in the 
hydrodynamical regime. However, the energy and the entropy slightly evolve 
in time in order to account for particles evaporation. 

Studying the evolution of different such systems with various initial masses, 
they have shown that, choosing a soft equation of state, the caloric curve may 
be explained as the temperature of the emitting source at the turning points 
of the collective expansion associated with various initial masses and exci- 
tation energies. They have also found that using a stiff equation of states 
reduces the amplitude of the considered monopole vibration leading to dif- 
ferent conclusions. 

The aim of our paper, is to examine these conclusions within microscopic 
model known to describe collective dynamics as well as single particle be- 
haviors in a quantum framework: the TDHF approach. In particular, we 
will test the validity of the breathing mode scenario by comparison with the 
complete TDHF dynamics. 



3 TDHF simulations 
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3.1 Mean- field approximation 



We consider the evolution of hot diluted spherical nuclei in the framework of 
mean-field theory. The simulation of the desexcitations of excited compressed 
spherical nuclei is already reported in the literature for nuclei initialized 
using a Constrained Hartree-Fock method (CHF). We use a scaling assump- 
tions in order to prepare the initial nucleus because this method allows to 
study a broad range of initial dilutions. In the cases where both methods are 
tractable, we have controlled that it gives results comparable with those of 



CHF [15 



We consider systems saturated in spin and isospin. Nucleons are moving 
independently in an average spherically-symmetric potential. We use the 
following parameterization for the effective mean-field potential 

U{p] = ~t p+ t3 p° +1 + c V 2 p+V c . (2) 

where the surface term, c V 2 p, can be related to the usual Skyrme parameters 
by 

5t 2 — 3£i 

For sake of simplicity we have taken forces with effective masses equal to the 
bare one. In order to reduce the numerical instabilities we have used the 
standard method which consists in replacing the surface term of potential 
(fj) by a folding product with a finite range function [[U| [L7j 



U [ P ] = \t'oP + ^4 P CT+1 + V Y ® p + V C (4) 



where Y is a Yukawa folding function 



Y{?) = (5) 



at the lowest order in the range a this expression is equivalent to the usual 
Skyrme functional provided that 

16 

t'o = t ~ y vra 3 V (6) 
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and that coefficients in the Yukawa folding functions are related by 



= (5t 2 - 3ti) = 
64vra 5 4vra 5 1 ' 

Finally, the direct Coulomb potential Vq is introduced in an approxi- 
mative way by giving an effective charge of +^ to each nucleon. Using 
this potential and a spherically symmetric density, each single-particle wave 
function can be separated into its radial, angular and spin-isospin part 

$ Q (f, a, r) = ^^Y lm (9, v)Xs(*)Xt(T) (8) 

where a represents all quantum number a = (n, I, m, s, t) , in which n is the 
energy principal number, (l,m) are the usual angular-momentum quantum 
numbers and (s,t) is the quantization of spin and isospin. If we consider 
initially a hot system, at temperature T, occupation numbers of various 
orbitals a are given by a Fermi-Dirac distribution 

1 

n n i = 7- — -x (9) 

exp + 1 

where the chemical potential is computed in order to get the correct number 
of particles. To be able to define particle orbitals at high temperature, a 
small external field Ar 2 has been added to the potential (A = 0.25 MeV/fm 2 ). 
Finally, the density of our spherical nucleus takes the particular form 

p(r,t) = 4j:(2l + l)n J K ^f , (10) 

n,l 

Before getting into details of the simulations, we will give some more 
details about these forces and about the associated Equation Of State (EOS). 

3.2 Forces and infinite nuclear matter properties 

Following [Q, we consider two parametrizations of the potential @ equiva- 
lent to the SI 1 1 and the SkM* forces (||18||). These two forces correspond 
respectively to a HARD and a SOFT EOS. 

Starting with the mean- field potential (|2]), the total energy per nucleon 
in the infinite, uniform, spin and isospin saturated nuclear matter is 

A A 8 16 v ; 
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where we have introduced the kinetic energy per nucleons for a cold gas 
of fermions this energy can be easily related to the density by 



Ek _ 
A ~ 

Therefore, the saturation density, po, is solution of 




(12) 



( d V\ n 2E ^ 3 + , + {*+!) nq x 

P*[f) =0 = 3T + 8^ + ^ (13) 

\ r / p=po 

which is nothing but the annulation of the pressure at the saturation point 
(at zero temperature). The infinite nuclear matter incompressibility is 
then given by 

\ r / p=po 

These three relations (|ll|), (0) and (14) provide an univocal relation between 



the saturation point properties p , -| and and the force parameters to, ^3 
and a. The parameter c does not influence the infinite uniform medium 
properties. 

If we compute the EOS associated with the parametrization ([|) we get 
similar expressions provided that we use the relation (||): t = t' Q + y7ra 3 Vo. 
The parameter Vq and a are related through the relation (0) but are not 
directly constrained by the EOS. However, for a finite system the relation 
between (Q) and (f|) is only valid at the lowest order in a. Therefore, this 
parameter should remains small. Following ref. [XT] we have fixed a to 
0.45979/m. Parameters values leading to same properties as SI 1 1 and SkM* 
for our parameterization are shown in table [I]. 

As an illustration, we have plotted the various isentropic EOS in infinite 
nuclear matter in figure [I] (Top). We have also plotted the isentropic and 
the isothermal spinodal regions defined as the region where the derivative 
of the pressure versus the density at constant entropy or temperature is 
negative. Therefore, in these regions the matter is mechanically unstable 
against density fluctuations: these are the so-called spinodal instabilities 
which might be responsible for the fragmentation of the system. 

Let us now come to the simulation of finite systems. These numerical 
simulations consist in two distinct parts: the first one is the initialization of 
abraded nuclei and the second one the TDHF dynamics. 



7 



3.3 Initialization of finite systems 

The first step is to initialize our system according to the equation (|). The 
method used in order to generate a hot and diluted initial source proceeds 
as follow: 

• First we look for the ground state of the nucleus with A4 = A — AAi 
nucleons. To do so we use the imaginary time method with a small 
constraint Ar 2 added to the mean field in order to get occupied states, 
\&h) ( h for hole states), and the associated one-body density p = 
J2h \®h) n h (®h\ where the occupation number are given by equation 
© at zero temperature. Diagonalizing the Schrodinger equation 

+ ^r4^ + U( P (r)) + Ar 2 ) R nl (r) = e nl R nl (r) (15) 
[ 2m dr z 2mr z J 

we can define the single particle states (nl). The small constraint Ar 2 
insures a limited number of orbitals for treating the continuum. In anal- 
ogy with the nuclear matter, we define the saturation density, p 
as the central averaged density in a sphere of 2fm. 

• Then, single particle states are occupied according to a Fermi-Dirac 
statistic ([D using a temperature T as a free parameter. This tempera- 
ture is defined iteratively in such a way that after the scaling described 
below, we get the excitation energy defined by Eq. ([!])[]. The chemical 
potential is determined requiring that the total number of particle is 
A — AAi. The initial entropy of the system is defined by 

S = -4j2(U+ 1) (n nl \ogn nl + (1 - n nl ) log(l - n nl )) (17) 

nl 

^or each AAj, we can define the excitation energy by subtracting the ground state 
energy of a cold nucleus of size A — AAi to the energy of the excited nucleus after rescaling. 
For all these calculations, we have used the expression of the energy corresponding to the 
potential (g) 

E = E kin + J /|t p(f)2 + it 3 p(f)('+a) + ^ J Y(r-r')p(r)p(r')dr'\df+E c (16) 
In this expression, E^in and Eq are respectively the kinetic and Coulomb energy. 
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• Finally we apply a scaling to each radial wave functions 

RT{r)=MR nl (-) (18) 

where M is a normalization factor. Then we compute the average 
density p (AA at the center of the nucleus (within a sphere of 2 fm) and 
we fix 5 in order to get the desired dilatation factor 77 = p ^f?\ = 0.8. 

After the initialization step we are following an ensemble of diluted nuclei 
with mass A — AAi and corresponding excitation energy E*. Each one is 
associated with an initial entropy A sample of various initial conditions 
are given m table |. Fi gure g presents various density profiles. On this graph 
one can spot only minor differences between the two forces. 

In order to get a deeper insight in the properties of finite systems, we 
have plotted the energy of the considered nuclei as a function of their central 
density: E = f(rf) (fig.fll]) middle and bottom). Each curve is obtained by 
scaling the wave-functions at fixed occupation numbers (which means that 
along these curves S = cie0). Middle part of fig. (JT|) corresponds to a 7g 7 Au 
for various entropies whereas bottom corresponds to masses and entropies 
reported in table |[ 

These curves can be compared to the infinite medium EOS shown in figure 
|] (Top). As far as the dynamics is concerned, these isentropic curves have 
only a meaning if the time dependent density profiles could be obtained one 
from another by a simple scaling in r-space. This is equivalent to consider 
the breathing mode as dominant and to neglect other fluctuations in density. 

Figure 1 illustrates the fact that finite systems have surface and Coulomb 
contributions to their energy. This modifies the saturation energy to 8-9 MeV 
(depending upon the mass of the nucleus considered) instead of 16 MeV found 
in the case of a neutral infinite system. Moreover, it changes the curvature of 
the isentropes i.e. it changes the effective compressibility modulus. However, 
this compressibility remains the main difference between the two ensembles 
of calculations associated with the two different forces. These isentropes will 
help us in the understanding of the dynamical evolutions since, in a naive 
picture, the collective vibration should oscillate around the minimum of the 

2 These relations between E , 77 and S can be considered as a finite system equation of 
states. 
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energy for the considered entropy and the turning points of the breathing 
mode should be simply obtained by requiring energy conservation 0. 

3.4 TDHF evolution 

Having defined the initial conditions as a given nucleus of size A t = A — AA^, 
a given excitation energy and the dilatation coefficient r] = 0.8, we let the 
system evolve through the TDHF equations f] 

^dRni(r,t) i-h 2 d h 2 l(l + l) rrrw .1 . . 

The potential term U [p] is the same as in equation (|J). The TDHF 
equations have the particularity to conserve occupation numbers (n n i) of 
single particle levels. 

|n- = (20) 

The entropy of the whole system (|T7| ) is thus conserved during the evolution. 
On the other hand the total energy of the system is also conserved. However, 
during the evolution, nuclei are evaporating particles and a tail at large 
distance will develop for the various wave functions R n i. We can try to 
investigate this process using the methods developed in which consist in 
splitting the r-space into two parts : 

• a sphere containing the initial nucleus, here we have considered a sphere 
of radius 15 fm. This part of the system will be called the nucleus in 
the rest of this article. 

• the rest of the space which will be considered as evaporated particles^. 

Therefore, due to the particle evaporation, the nucleus is continuously loosing 
mass, energy and entropy f\ We are thus in the same framework as in ref. @, 

3 Note that, picture ([!]) depends upon the size of the sphere we have taken in order to 
define the density ( here r=2 fm). This will be discussed in the following. 

4 For the evolution, we have taken a step of time At = 0.75/m/c. We have discretize 
the r-space in step of size Ar = 0.2 fm, the total size of the r-space being 300 fm, this size 
is big enough to avoid bouncing of evaporated particles against boundaries. 

5 We are considering very large space in order to avoid the reflections of the wave packets 
on the boundaries of the considered r-space. 

6 In general at very low excitation energy the nucleus can even gain entropy because 
of the dissipation of the collective motion. 
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T§| but with a microscopic model. However, the TDHF approach is neglecting 



the effect of the collisions and this might induce differences between our 



simulations and those of ref. 0, These differences are expected to be not 



too important because in a self-similar expansion of infinite nuclear matter, 
the gain and the loss term due to collisions cancel, so that the mean field 
dominates the dynamics. In a finite system, the collisions do not exactly 
cancel as it is known from the calculation of the small amplitude motion ]20| . 
This should be kept in mind when the various models are compared. 



4 Discussion of the expansion dynamics 

4.1 Evolution of global quantities: evaporation and 
monopole expansion 

In this section, we study the evolution of the various considered nuclei. In 
fig. (|3p, we have plotted the evolution of masses and entropies of the systems 
reported in table |[ All these quantities are computed within a sphere of 15 
fm. From this picture we can see that the smaller are the initial masses, the 
larger are the energies and entropies and so are the number of evaporated 
particles. Indeed, one can observe that the evaporation is stronger for light 
highly excited systems with a rapid decrease of the mass and the entropy. 
Moreover, in TDHF approach the evaporation is treated in quantum me- 
chanics, particle may be reflected several time by the potential well before 
being ejected from the nucleus. Therefore, the particle evaporation appears 
rather slow. 

This monopole motion can be clearly seen in figure [| which displays the 
time evolution of the central densities for different initial conditions. From 
this figure, as expected, collective motion is slower and presents a larger 
amplitude motion in the case of a soft EOS. 

4.2 Discussion of the Collective dynamics 

In these figures, we do not observe the fast expansion towards low density 
regions. Indeed, at moderate energies (T < 5 MeV), we observe that the first 
minimum (i.e. the first turning point) appears to be close from the initial 
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dilution (see for example A = 191). At low excitation energies, the small 
dilution of the abrasion-ablation stage is followed by a recompression phase. 
Then, when the system expands again, the average density may reach values 
lower than the initial one. 

Conversely, in a simple description assuming a collective dynamics at 
almost constant entropy and mass, this behavior should be normally under- 
stood considering the energetics of isentropes as shown in figure [|. Indeed, 
schematically, if the system is initially with an excitation energy, Eq, and 
an entropy, So, it will expand or contract till it has reached the other point 
of the isentrope S = Sq which cut a horizontal line corresponding to the 
energy E = E . This schematic procedure is illustrated in figure |]. Values 
of turning points obtained with this simple construction as well as the one 
of ref [U are compared to those extracted from the simulation in table [|. 

For A = 150 and A = 111, one can see only little differences between 
results of the schematic construction and those of the exact dynamics. How- 
ever, they are different for A = 191 and A = 89. Considering table ^ and 
figure |], several comments should be made: 

• i) For high excitation energies (small masses), the schematic construc- 
tion based on energy and entropy conservation predicts a total vapor- 
ization of the nucleus whereas in the dynamical simulations a small 
residue survives. This discrepancy can be accounted for by considering 
the energy and entropy variation due to the evaporation of particles 
(see fig.|). 

• ii) for low excitation energy (for example A = 191), evaporation is 
negligeable (S=cte, T=cte (fig. ^)) and the result could be directly 
compared to the schematic construction. However, we observe that 
TDHF predictions differ from schematic predictions. In the schematic 
model, from figure |6|, starting with a nucleus of mass A = 191 at 
f] = 0.8, it appears impossible to reach a dilution lower than the initial 
one since the diluted nucleus is predicted to undergo a compression 
right after the beginning of the calculation. In TDHF simulations, the 
system is able to reach lower densities than the initial one by creating 
a hole at the center of the nucleus. Moreover, instead of undergoing 
damped monopolar oscillations, the nucleus reaches a much lower value 
of central density at the second dilatationF]. 

7 Note that, this effect is not due to the Coulomb interaction and so is not analogous 
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As far as the comparison of TDHF simulation and the results of 0] is 
concerned, the values of the turning points reported in table f| appears to 
be different. At high initial temperature (T > 5 MeV) and for masses lower 
than 150, the amplitude of the expansion is larger in the collective model 
than in TDHF. At low temperature, the TDHF expansion is more complex 
than a breathing mode picture. In this picture, the central density is taken 
to be the only relevant variable since a self similar expansion is explicitly 
assumed. In fig. |5|. we show a typical example of density profile for different 
times of the evolution for A = 191. From this picture, it is clear that hollow 
structures are created during the evolution and that one density profile at 
a given time could not be obtained from another one by a simple scaling. 
The existence of exotic shapes which are created during the expansion, was 
already suggested in several works ||. The main novelty found here as we 
will see, is that these structures are consequences of the beating of intrinsic 
collective modes[[[3|]. In order to get a deeper insight, we will focus on the 



collective dynamics in the following^. 
4.3 Collective dynamics 



In this section, we have generalized the analysis performed in ref . |L3[ to many 
different initial conditions. In TDHF calculation, no assumption is made on 
which collective degrees of freedom could develop during the evolution. 

In order to better understand the complex dynamics of the density p(r, t), 
we have plotted in the fig.|7| its Fourier transform p(r, uj) performed over 1500 
fm/c. This corresponds to the expected breathing vibration and indeed the 
radial dependence of the density variation p(r, u) is well fitted by the usual 
Tassie transition density [|12]1 . 

In fig.0, we observe several other collective higher at higher energies. For 
7] = 0.8, we have many additional waves at frequencies 17, 20, 25 and around 
40MeV. The collective motion located at 2 * 12.5 ~ 25MeV corresponds to 
the two-phonon excitation of the breathing mode vibration. The other peaks 
at 17, 20 and around 40MeV could be associated to other collective degrees 
of freedom. 



to the phenomenon reported by Borderie et al in ref. |9[, since removing the Coulomb 
potential from our simulations do not modify the observed behavior. 

8 In the rest of the paper, we will discuss only the soft-EOS since the effect we are 
discussing here are generic and, from the qualitative point of view, do not depend upon 
the compressibility of the force used in the simulation (see fig. [|). 
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In order to study either a possible anharmonicity or a coupling between 
modes, we have investigated the role of the vibration amplitude changing 
the initial conditions. In a pure harmonic picture, frequencies should remain 
constant and amplitudes of responses should vary linearly with the initial 
perturbation. Focussing on the middle and bottom part of fig. [7|, we see first 
that the frequencies are almost constant but the amplitude of the response 
do not depend linearly on the initial perturbation. Therefore, the nuclei 
presents non-linearities and mode couplings but they are not strong enough 
to destroy the collective motions. This means that we are still far from a 
chaotic regime. 

In this chapter, we have pointed out that not only the breathing mode 
is excited in our TDHF calculation but also other modes. These modes are 
highly non-local in r. This induces a complex dynamics and in particular 
may lead to the creation of hollow structure. Up to now, many different 
macroscopic models || |], |Tl| have been used in order to extract information 
about the break-up of highly excited nuclear systems. In these approaches, 
the breathing vibration is taken to be dominant and the expansion is assumed 
to be a time-dependent self-similar expansion of the nucleus. However, this 
picture fails to reproduce our calculation which indicate that we must extend 
our collective phase space to include more degrees of freedom. In the next 
chapter, we will address this question and illustrate a possible generalization 
the EOS of a finite nucleus in order to include hollow structures. 



5 Inclusion of Hollow shapes of the nucleus. 

During the evolution, we have seen that the nucleus could reach low density 
region by creating a hole at its center. In order to study this type of shape, 
we have parameterized the density profile as follow: 

p(r t) = A ( pi ■ a ^ 3 + eilM. ] , , i ) 

P[,) \l + e W (^^) + l + exp^M^' 



a-2 



In this expression, A is fixed in order to conserve the number of particles. 
The coefficients pi, Ri and ai are fixed by fitting the density obtained during 
the evolution of the nucleus A = 191 at time t = 108/m/c (see figure [5]). 
Pi, Ri and a\ are respectively equal to 0.0826/m~ 3 , 3.09/m and 0.598/m 
whereas p 2 , Ri and a 2 are equal to 0.1354/m -3 , 6.98/m and 0.528/m. With 
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this parameterization, the information contained in the density profile is now 
reduced to only two coefficients a(t) and (3{t). If a(t) is less than zero, the 
nucleus presents a hole at its center, otherwise, it presents a bump. Since 
P2 is much bigger than pi, for low values of a, the parameter (3 could be 
approximatively seen as a global scaling of the nucleus (/3 > 1 corresponds 
to a dilatation whereas (3 < 1 is a compression). 

In figure [5], we have traced results of the best fit we have obtained for the 
different density displayed. Values of a and (3 are reported in table ||. 

We can see that this parameterization seems to be particularly suitable 
and represents well the different density profiles, so that the time dependent 
evolution could be accurately replaced by the evolution of the two parameters 
a and (3. Since we now include a more general class of density shapes, this 
phase-space is a generalized approach of the breathing mode picture, the two 
variables being now interpreted as two collective coordinates. 



5.1 Generalization of the zero temperature equation 
of state. 



Using the parameterization fl2lD we can define the EOS in the (a, (3) phase 



space by computing the energy at zero temperature for different a and (3 



parameters. Since eq. (|2l|) defines only the density profile and not the wave- 



functions, we have used the Extended Thomas Fermi formalism described in 



181 where two additional terms are added to the usual Thomas Fermi kinetic 



E K (r) = E TF (r) + ( ^^y^ + ^Ap(r) ) (22) 





energy 



In this expression, Etf take the form (|T2|). The total energy is calculated by 
incorporating the potential (§) and integrating over the r-space. 

By varying the parameters, we can have the generalized EOS: E = 
E(a,j3). This EOS is represented in fig.|8]. The graphic of a self-similar 
picture displayed in figure |^ for A = 191 could be recovered by cutting the 
graphic with a surface a = f{(3) (which could be roughly approximated by 
a ~ cte. In this picture, we can see different potential wells with different 
minima with almost degenerated energies separated by small barriers (of the 
order of 1 — 2 MeV). These three minima correspond respectively to a nucleus 
with a hole (a < 0), a standard Wood-Saxon shape (a = 0) and a nucleus 
with a bump at the center (a > 0). 
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In order to study the evolution of the system in this E(a, (3) phase-space, 
we have performed a fit of the density profile at each time. Results is shown 
up to 108/m/ c in figure ^. The path in this phase-space appears complex and 
cannot be reduced to a scaling of the (3 variable. This is another indication 
of the fact that we cannot reduce the TDHF dynamics to a breathing mode 
picture. During the evolution, we see that the nucleus could go in many 
different configurations in this phase space. 

In summary, the dynamical expansion as the one studied by 0] assumes 
a predominance of the breathing mode. In a microscopic calculation, other 
collective vibrations can be excited during the dynamical evolution. This 
lead to a complex dynamics for hot and compressed (or dilated) systems 
for moderate energies. As a direct consequence, the dynamics should be 
discussed in an enlarged phase-space of several collective degrees-of-freedom. 
The illustrative pictures and |^) seem to indicate that energetically hollow 
and uniform nuclei should be both taken into account on the same level. Note 
however, that for initial temperature greater than 5 MeV, many orbitals are 
occupied and levels crossing seems less important and we do not observe this 
effect anymore. 

Finally, we want to mention that the theoretical determination of an ef- 
fective temperature during the evolution of the expanding source, is often a 
hard task and is often not possible since the system is only equilibrated at the 
initial time. During the dynamics, the nucleus is strongly out-of equilibrium 
and no energy and occupations numbers of levels could be defined on the 
same time (no approximative Fermi-Dirac distribution could be defined f§. 
However, we can imagine another way of defining the temperature. Indeed, 
it is possible to follow the density and the energy of the excited system with 
time. In a breathing mode picture, we have a univocal correpondance be- 
tween these two parameters and the temperature: i.e. the EOS E = E(p, T). 
However, this picture breaks in our calculation since we have to consider at 
least E = E(a,/3,T) and in fact two parameters are often not enough. In 
conclusion, the determination of an effective temperature in models as the 
one we are considering, is not well-defined as far as no equilibration process 
is added on top of the mean-field. 



9 Indeed, during the evolution the Hamiltonian h and the density p do not commute. 
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6 Discussion 

We have used a quantum microscopic simulation in order to discuss the inter- 
pretation of the caloric curve recently proposed by Papp and N6renbergf|. 
In their work the experimentally observed caloric curve may be interpreted as 
a line of turning-point reached in the low density region during the monopole 
expansion of excited nuclei where the nucleus could break. Using analogous 
initial conditions, we have pointed out that many differences between a mi- 
croscopic approach and the collective model. Indeed at low and moderate 
energies, the TDHF dynamics indicates that different collective vibrations 
could develop on top of the self-similar vibration. This collective behavior 
may lead to the formation of hollow structures. Since the density come from 
a superposition of many modes, we see that the lowest density region could 
be reached after a longer time (see fig. |j) than the first breathing of the 
nucleus. It should be noticed that experimental DATA of photon detection 
seems to indicate such a behavior ||21||. The presence of many collective modes 
in the dynamic has led us to generalize the simple breathing mode picture. 
By a suitable parameterization of the density, we have defined a generalized 
Equation Of State. This new approach, include the possible existence of 
exotic shape in the dynamic. As we have seen, the nucleus has a complex 
path in this new phase-space. On the other hand, at high initial excitation 
energy, when the breathing collective mode dominates, the amplitude of the 
monopole oscillation is also different from the one reported in[Q. Since our 
simulations demonstrate that the collective dynamics of an abraded nucleus 
could be complex and depends on the model used. The possible interpre- 
tation of the caloric curve as a line of turning-point has to be considered 
carefully. 
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Soft EOS (SkM*) 


Hard EOS (SIII) 


f (MeV) 


-15.77 

1 


-15.85 
1 


Po (fm- 3 ) 

to (MeV fm 3 ) 

t 3 (MeV fm 3 ( 1 +' T )) 

a (fm) 

Vb (MeV) 

i^oo (MeV) 


6 

0.16 
-2191.73 
18818.8 
0.45979 
-461.07 
198.88 


0.145 
-435.41 
17258.8 
0.45979 
-401.77 

368.4 



Table 1: Parameters and properties of the two interactions used in this article. 



A 


Si Soft EOS (SkM*) 


St Hard EOS (SIII) 


191 


0.44 


0.37 


187 


0.38 


0.35 


183 


0.45 


0.51 


177 


0.63 


0.51 


172 


0.79 


0.71 


165 


0.97 


0.92 


158 


1.14 


1.12 


150 


1.33 


1.31 


142 


1.52 


1.50 


134 


1.59 


1.55 


120 


1.80 


1.82 


111 


2.0 


2.06 


98 


2.46 


2.49 


89 


2.79 


2.80 



Table 2: Examples of initial conditions considered. For each initial masses 
and for the two considered forces, initial entropies are given. 
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SOFT EOS 



HARD EOS 




Figure 1: Relations between the energy per nucleon and the central density. 
Top parts correspond to the infinit medium EOS, each curve is drawn at 
constant entropy and the entropy variation between two adjacent curves is 
AS/ks = 1. In such a case both the isothermal (dashed) and the isentropic 
(solid) spinodal are displayed as thick lines. Middle part corresponds to a 
79 7 Au nucleus while the bottom part is associated with the various systems 
considered for the evolution. In the bottom part, curves correspond to de- 
creasing masses (equivalently to increasing excitation energy) from bottom 
to top of the graphic. More detail is given in the text for the determination 
of this picture in finite nuclei. For all graphs, left and right part correspond 
respectively to a soft and a hard EOS. 



A 

Ai 


Po [Ai) Soft EOS (5km ) 


Po [Ai) Hard EOS {bill) 


191 


0.142 


0.139 


187 


0.142 


0.139 


183 


0.143 


0.140 


177 


0.147 


0.142 


172 


0.150 


0.144 


165 


0.155 


0.147 


158 


0.160 


0.150 


150 


0.165 


0.153 


142 


0.170 


0.157 


134 


0.155 


0.145 


120 


0.143 


0.140 


111 


0.143 


0.140 


98 


0.154 


0.148 


89 


0.163 


0.154 



Table 3: Saturation densities for the various initial masses. 



Soft EOS {SkM* 



Hard EOS (SIII) 




Pmin/ PO 


Pmin/ PO 


Pmin/ PO 




Pmin/ PO 


Pmin/ PO 


Pmin/ PO 


from [4] 


schem 


exact 




from [4] 


schem 


exact 


0.68 


1.45 


0.408 




0.75 


1.32 


0.67 


0.4 


0.82 


0.764 




0.52 


0.86 


0.80 


0.0 


0.6 


0.77 




0.23 


0.89 


0.786 


0.0 


0.0 


0.25 




0.0 


0.0 


0.494 



Table 4: Comparison between the turning point observed dynamically and 
those predicted using a schematic self-similar isentropic dynamics at constant 
energy (see figure 6). Turning point given in [4] for the same forces are also 
reported. 



t (fm/c) 


a(t) 


P(t) 





0.60 


1.07 


18 


-0.80 


0.98 


45 


0.76 


0.89 


108 


-1.00 


1.00 



Table 5: Different values of a(t) and f3(t) obtained by fitting the density 
profile at divers time of the expansion of the nucleus A = 191. 
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HARD 




0.00 



0.0 2.5 5.0 7.5 10.0 12.5 0.0 2.5 5.0 7.5 10.0 12.5 15.0 

r(fm) r(fm) 



Figure 2: Typical examples of initial density profiles for various initial masses: 
Left (soft EOS) and right (Hard EOS). 
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Figure 3: Typical examples of time evolution of the mass, and entropy for 
the various initial masses considered. 
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Figure 4: Typical examples of the time evolution of the central density for 
various initial masses (A = 191, 150, 111 and 89.4). 
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0.15 — 
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0.20 



0.10 



0.05 
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0.00 
0.20 — 
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0.05 — 



-H- 
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t=108 fm/c 




10 
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Figure 5: Density profiles taken when the central density reaches its extrema. 
Figures corresponds to an initial mass A = 191 and time t — 0, 18,45 and 
108 fm/c from top to bottom, which corresponds to turning points of the 
central density. Solid line corresponds to density profiles obtained in TDHF 
calculation using a soft EOS. Dashed lines corresponds to fit obtained with 
two wood-saxon (see text). 
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SOFT EOS 




Figure 6: Shematic construction of the collective turning point using the 
isentropic curve in the excitation energy versus density plot, asking for both 
a conservation of the excitation energy and the entropy. Initial conditions 
are reported as crosses whereas first turning points are denoted by circles. 
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Figure 7: Fourier transform p(r,u) (Left) p(r,t). From top to bottom three 
different initial compression/dilatation factor are considered: rj = 0.8,1.0 
and 1.2 for a nucleus of mass A = 191. 
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Figure 8: Representation of a generalised EOS which include hollow struc- 
tures. The energy is reported as a function of a and (3. a greater than 
corresponds to nucleus with a bump whereas a less than one corresponds to 
a nucleus with a hole at its center. On the other side, f3 > 1 is equivalent to 
a compression and (3 < 1 means a dilatation. 



29 




1.0 -0.5 0.0 0.5 1.0 1.5 



Figure 9: Dynamical Path of TDHF simulation in the phase-space. 
Contour plot (dashed lines) of the energy are also displayed (contour corre- 
sponds to -8, -7, -6, -5 and -AMeV . 
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